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Attosecond XUV vortices from high-order harmonic generation 

C. Hernandez-Garci'a/’[^ A. Picon,^ J. San Roman,^ and L. Plaja^ 

^ Grupo de Investigacion en Optica Extrema, Universidad de Salamanca, E-37008, Salamanca, Spain 
^Argonne National Laboratory, Argonne, Illinois 60439, USA 
(Dated: July 3, 2015) 

We present a theoretical study of high-order harmonic generation (HHG) and propagation driven 
by an infrared field carrying orbital angular momentum (0AM). Our calculations unveil the fol¬ 
lowing relevant phenomena: extreme-ultraviolet harmonic vortices are generated and survive to the 
propagation effects, vortices transport high-OAM multiple of the corresponding 0AM of the driving 
field and, finally, the different harmonic vortices are emitted with similar divergence. We also show 
the possibility of combining 0AM and HHG phase-locking to produce attosecond pulses with helical 
pulse structure. 
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Helical phased beams, also called optical vortices, are 
structures of the electromagnetic field with a spiral phase 
ramp about a point-phase singularity. The phase wind 
imprints an orbital angular momentum (0AM) to the 
beam, in addition to the intrinsic angular momentum 
associated with the polarization m\- As light-matter 
interaction is inherently connected with the exchange of 
momentum, 0AM can also be transferred to atoms ms], 
molecules [3 IB], or an ensemble of atoms . Opti¬ 
cal vortices have potential technological applications in 
optical communication dlH], micromanipulation [IB], 
phase-contrast microscopy (TT] [TH|, and others [T^. In 
the recent years an interest has burgeoned in imprint¬ 
ing phase singularities in the XUV/X-ray light generated 
both in synchrotron and X-ray free-electron laser (XFEL) 
facilities pOl - E^ . In this short-wavelength regime one can 
drastically reduce the diffraction limit, as well as exploit 
the selectively site-specific excitation, with an important 
impact in microscopy [221 [22] and spectroscopy |2S] . 

Phase singularities can also be imprinted to shorter- 
wavelength light using high-order harmonic generation 
(HHG), as it was recently demonstrated by Ziirch and 
coworkers producing vortices in the extreme ultraviolet 
(XUV) [25]. The confluence of 0AM and HHG con¬ 
stitutes an extraordinary promising perspective. For 
instance, the phase twist is not imprinted directly to 
the short-wavelength radiation, but to the fundamental 
field. Therefore, it requires a single setup (diffractive 
mask, for instance) to imprint 0AM to the fundamental 
field, that will be subsequently transferred to the rain¬ 
bow of harmonic wavelengths by non-linear conversion. 
On the other hand, high-order harmonics have extraor¬ 
dinary temporal coherence qualities 1211123 , a unique 
feature that has triggered a revolutionary metrology tool 
for the temporal characterization of ultrafast processes at 
the atomic scale (both spatially and temporally) [221EB] ■ 
Two decades ago, it was demonstrated that the higher 
part of the HHG spectra can be used to synthesize short 
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XUV pulses of attosecond duration [2IH32]- It is, there¬ 
fore, an appealing possibility if such attosecond pulses 
can be synthesized with 0AM. In addition, the new HHG 
generation schemes, based on the present development of 
mid-infrared (mid-IR) laser sources, show that keV X- 
ray radiation can be obtained from the HHG process in 
a table-top system [55], as well as that the zeptosecond 
era could be closer than expected [5S] . Therefore, there is 
no fundamental obstacle for up-shifting the XUV 0AM 
beams to extremely short temporal structures, or to gen¬ 
erate them in the soft X-ray regime. 

A surprising finding of the recent HHG-OAM experi¬ 
ment [26] is that the topological charge of the harmonic 
vortices is nearly equal to one, i.e. that of the funda¬ 
mental field. This is counterintuitive in terms of the 
present understanding of HHG with intense fields, in 
which the phase of the harmonics scales roughly with 
the harmonic order. At present, it is recognized that 
a detailed theoretical treatment, including propagation 
effects, is needed to elucidate some of the fundamental 
aspects of the HHG 0AM conversion m- For instance, 
to confirm the increase of the topological charge with the 
harmonic order and, thus, to corroborate that the origin 
of the lower topological charge of the harmonics detected 
in [26] should be attributed to parametric instabilities of 
the nonlinear propagation and not to the fundamental 
HHG process. It is presently uncertain that high-charge 
vortices might be obtained, even if no parametric insta¬ 
bilities are present, as they would have to survive phase¬ 
matching effects during propagation. Finally, it is also 
an open question whether attosecond HHG 0AM pulses 
can be generated and be, also, resilient to propagation. 

In this Letter we present a pioneering theoretical study 
of the HHG process from an IR field carrying 0AM of 
topological charge £, in the low intensity regime. In this 
regime, the nonlinear propagation instabilities, which 
could affect drastically the observations [26] j can be ne¬ 
glected. We prove that: 1) XUV harmonics vortices are 
generated and survive to the propagation effects, 2) each 
harmonic has a topological charge of q£, q being the har¬ 
monic order, 3) all the harmonics are emitted with similar 
divergence and, 4) attosecond pulses carrying 0AM can 
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be generated and also survive propagation. 

We compute harmonic propagation using the electro¬ 
magnetic field propagator [35] ■ The HHG emission of the 
single-atom sources is computed using the SFA-I- method, 
an extension of the standard strong field approach with 
good quantitative accuracy [39|. Every point of the tar¬ 
get (gas cell or gas jet) is treated as source of an elemen¬ 
tary wave generated by an ensemble of atoms, which is 
then propagated to the detector. The final field at the 
detector is, thus, the coherent addition of these elemen¬ 
tary waves. One of the advantages of this method is that 
is well-fitted to compute high-order harmonic propaga¬ 
tion in non-symmetric geometries, therefore, it is spe¬ 
cially suited for computing HHG driven by beams car¬ 
rying 0AM. The method has been successfully used for 
describing regular HHG with near- and mid-IR lasers, in 
good agreement with experiments |36l [40] . More infor¬ 
mation on this method and the SFA-I- may be found in 
Refs. [311133] . 

In order to implement the 0AM beams, we consider 
the Laguerre-Gaussian (LG) modes in the paraxial ap¬ 
proximation [T], i.e. the beam propagates in a well- 
defined direction (in our case the z-axis) as a plane wave 
modulated by the slowly-varying transverse amplitude 
(written in cylindrical coordinates) 
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W{z) = Wo^/l + {zjzoY is the beam width, where Wq 
is the beam waist, Wq = \f\zol^ (-^o being the Rayleigh 
length); ELiz) is the wavefront radius of curvature, given 
by R{z) = z[l-|-(zo/z)^]; ((z) is the Gouy phase, which is 
given by C(z) = -{\^\ -I-2p-|- 1) tan“^(z/zo); and, Lp'[a;] 
are the associated Laguerre polynomials. The indices 
i = 0, ±1,±2,... and p = 0,1,2,... correspond to the 
topological charge and the number of radial nodes of the 
mode, respectively. 

We consider as a fundamental field a beam carrying 
0AM with .£ = 1, or £ = 2, and p=0, with a beam waist 
of Wo = 30 /rm, and hence a Rayleigh range zo=3.5 mm. 
The amplitude of the field, Eo is chosen to give a peak 
intensity at focus of 1.4 x 10^"^ W/cm2. The laser pulse 
is assumed to be a sin^ envelope of 5.8 cycles (15.4 fs) 
FWHM and 800 nm wavelength. In Figure [^a,b) we 
present the intensity and phase profiles of the LGi^ mode 
at the focus position for the above parameters. As it 
can be observed in plot (b), the term imprints an 
azimuthal phase variation on the beam from —tt to tt. 

In our simulations, the 0AM beam is focused into an 
argon gas jet, which is directed along the x-axis, and 
modeled by a Gaussian distribution along the y and z 
dimensions (whose FWHM is 500 ^m), and a constant 
profile along its axial dimension, a;, with a peak density 
of 10^^ atoms/cm^. The gas jet is located after the focus 
position, because in this region the Gouy phase and the 
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FIG. 1. (a) Intensity and (b) phase transversal profiles of the 
Z/Gi,o rnode of the fundamental beam at the focus position, 
with a beam waist of Wo = 30 pm. The amplitude Eo in Eq. 
Q, is chosen to give a peak intensity of 1.4 x W/cm^ at 
the focus. (c,e) Intensity and (d,f) phase angular profiles for 
the 17th-23rd harmonics. Note that the intensity profiles are 
given in logarithmic scale. The resulting topological charge 
can be obtained from (d) and (f) resulting in £ = 17 for the 
17th and £ = 23 for the 23rd harmonic. 


intrinsic phase of the harmonics tend to compensate each 
other, resulting in optimal longitudinal phase-matching 
conditions [43l|34]. In addition, as mentioned above, the 
pulse intensity and the gas density are chosen such that 
the nonlinear effects of the propagation can be neglected 
both for the IR as well as the XUV fields. 

In Fig. [^c-f) we present the intensity-phase angular 
profiles of the 17th (c,d) and 23rd (e,f) harmonics gener¬ 
ated with the LGifi mode. Two main conclusions arise 
from these plots. First, the radius of the annular inten¬ 
sity distribution is similar for the two selected harmonics, 
thus, they are emitted with similar divergence. In fact, 
similar divergence emission can be observed for the whole 
harmonic spectra -see Fig. ia)-. Secondly, the phase- 
plots show that the charge of the gth-ovdev harmonic is 
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g, which, as mentioned above, is expected from the HHG 
theory m- The harmonics generated with the TG 2 ,o 
mode (not shown) fulfill also these properties, appear¬ 
ing with a similar divergence angle and the gth-harmonic 
having 2q topological charge. 



Divergence angle (mrad) 


29 

27 

Fundamental; £ = 2 


25 

U 

a> 



"2 23 
0 

•a 


.y 21 
0 



1 19 

ro 

^ 17 




Divergence angle (mrad) 


FIG. 2. Angular distribution of the harmonic spectra gener¬ 
ated by an optical vortex of charge (a) £ = 1 (LGi,o) and (b) 
1 = 2 (LG 2 ,o). The parameters for plot (a) are the same as in 
Fig. 0 whereas for plot (b) we have changed the topological 
charge, 1 = 2, and the amplitude Eg in Eq. 0 has been 
increased to give a peak intensity at the focus of 1.4 x 10^"^ 
W/cm^. The radii of the generated XUV vortex is similar 
for all the harmonics, being longer for the higher topological 
charge. The spectra are presented assuming a A1 filter plate 
of 500 nm in thickness. 


Figure [^a) shows the angular profile of the spec¬ 
trum diverging along the y-axis (transverse to propaga¬ 
tion) for the same parameters as in Fig. 0 We can 
observe that the harmonics in the plateau region ex¬ 
hibit a double-fringe profile corresponding to the vor¬ 
tex intensity-distribution, whose radius is similar for all 
the harmonics, a feature that was also observed in [50]. 
For a fundamental beam of topological charge £ = 2, 


in particular a LG 2 ,o rnode, we also observe that the 
high-harmonics are emitted as XUV vortices with simi¬ 
lar radii -see Fig. ib)-, being the radii longer than 
those obtained with £ = 1. Therefore, HHG leads to a 
perfect vortex generation process in terms of its applica¬ 
bility |45| , as all these XUV vortices of topological charge 
q£ are emitted with similar size. 

Let us now look at how the XUV vortices are emit¬ 
ted temporally. One of the most exciting perspec¬ 
tives of high-order harmonic generation by intense lasers 
is the possibility of synthesizing XUV pulses of sub¬ 
femtosecond duration [3T|. An attosecond pulse train (1 
as = 10“^® s) is obtained by the selection of the higher 
frequency part of the harmonic spectrum that conform 
the plateau region |53|. For the correct synthesis, the 
spectrum should approximately satisfy the two following 
conditions: on one side its structure should approach to 
that of a frequency comb, in which the harmonic inten¬ 
sities are similar; on the other side, the relative phase 
between the harmonics should be nearly constant (phase 
locking) |3T|. These two conditions are well satisfied in 
the typical harmonic spectrum generated during the in¬ 
teraction of an intense field with an atom. Fortunately, 
HHG driven by 0AM beams leads to the generation of 
XUV 0AM beams with similar divergence, therefore al¬ 
lowing for the synthesis of an helical attosecond pulse 
train. 

In figure 0a) we present the temporal evolution of the 
high-harmonic signal shown in Figs. 0 and 0a). As 
the fundamental laser pulse consists in many cycles - 
5.8 cycles (15.4 fs)-, an helical attosecond pulse train is 
obtained, i.e. an attosecond pulse train delayed along 
the azimuthal coordinate according to the phase varia¬ 
tion of the fundamental LGi^ beam. Note that the time 
separation at a fixed angle between the pulses of the he¬ 
lical structure is half the period of the driving held (1.33 
fs). On the right side of Fig. 0a) we show the attosec¬ 
ond pulse train obtained at different azimuthal angles 4>. 
Note that although each high-harmonic 0AM beam is 
generated with £ = q, as we are selecting many harmonic 
orders (see Fig. 0, the pulse train structure remains 
similar. 

It is also possible to obtain an isolated attosecond pulse 
using few-cycle driving pulses |32j [34] , since high-order 
harmonics are then generated in a single rescattering 
event. For that purpose we show in hgure0b) the tem¬ 
poral structure of the XUV 0AM beam driven by a few- 
cycle pulse -1.4 cycles (3.8 fs)-. The 0AM imprints a 
different carrier-envelope phase (CEP) in the fundamen¬ 
tal pulse over the azimuthal angle 4>, and thus, the helical 
attosecond structure varies from an isolated attosecond 
pulse, generated around (f) = 60°, to a double pulse struc¬ 
ture, around (j) = 120°. As a consequence, we obtain in a 
single-shot the map of all the attosecond pulse structures 
for different CEP values. We believe that these helical 
attosecond structures, obtained in a single-shot, are a 
powerful tool for pump-probe experiments with attosec¬ 
ond resolution. 
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FIG. 3. Temporal evolution of the high-harmonic signal driven by (a) a multi-cycle pulse -5.8 cycles (15.4 fs) FWHM- and 
(b) a few-cycle pulse -1.4 cycles (3.8 fs) FWHM-. The XUV 0AM is emitted in the form of an helical attosecond structure, 
whose detailed form is presented on the right side for different azimuthal angles (f>'. 0°, 60°, 120° and 180°. 


In conclusion, we have presented theoretical calcula¬ 
tions of the HHG process driven by an intense IR beam 
carrying 0AM with topological charge £. We have ob¬ 
tained XUV harmonic vortices with topological charge 
of q£, q being the harmonic order. In addition, we have 
shown that all the harmonics are emitted with similar 
divergence and that these structures are robust under 
propagation. We have also proven that the highest order 
harmonics can be used to synthesize 0AM attosecond 
pulses, that exhibit an attosecond helical structure. We 
believe that the combination of the properties of optical 
vortices carrying 0AM with the spatio-temporal charac¬ 


teristics of high-order harmonic generation and propaga¬ 
tion opens a new perspective in ultrafast science. 
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